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Spacing interval between principal tree windbreaks 
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Abstract: Relative windspeed reduction was measured behind nine relatively narrow, homogenous tree windbreaks with poros¬ 
ities between 0.13-0.33, and behind 28 combinations of model stubble barriers representing 25 different optical porosities 
(0.00-0.80). The optimum porosities observed were 0.25 and 0.13 for tree windbreaks and stubble barriers respectively. Based 
on the relationship between windbreak structure (optical porosity) and wind reduction, the chief indices for determining spacing 
interval, i.e., the windbreak structure index (5) and the parameter of microclimate, represented by the problem wind (Lrp), were 
determined. Additionally, investigations on shelterbelt trees were carried out, and stem-analysis techniques were used, to de¬ 
velop a method for determining the mature height of tree windbreaks ( H 0 ). Optimal spacing intervals between windbreaks could 
be predicted from the indices of a given windbreak structure, percentage of reduction of windspeed desired and tree growth 
model. A hypothetical example for determining the spacing interval of principal poplar windbreaks is given at the end of this 
paper. The results can be applied not only to tree windbreak design but also to other plant materials and artificial barriers for 
wind protection. 
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Introduction 

The basic principles of wind control indicate that wind¬ 
breaks planted or placed in the path of air flow act as sur¬ 
face barriers, causing an upward diversion of the air current 
(Woodruff 1956; Skidmore and Hagen 1977; Maki 1985; 
Zhu et al. 1998, 2001), and this diversion is accompanied 
by a drag on the wind at certain heights of the windbreaks. 
These combined effects lessen the force of wind on the 
original ground surface, lower the prevailing surface veloc¬ 
ity (Zhu et al. 2000), and create an area of relatively calm 
airflow within the zone of influence of the windbreaks. This 
zone of influence is limited, extending leeward to a distance 
from several to more than 30 times the height of the wind¬ 
break (Woodruff 1956; Cao et al. 1981). However, the re¬ 
quirements for protection usually exceed this limited influ¬ 
ence. In the case of windbreaks, protection can be 
achieved by establishing a system of properly oriented 
windbreaks. The most common design parameter for a 
system of windbreaks is the spacing interval (Cao 1983). 

Many researchers (Woodruff 1956; Cao 1983; Tibke 
1988; Ticknor 1988; Finch 1988) use the following simple 
equation to determine the spacing interval. 
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SlP = cH 0 (1) 

where SIP is the spacing interval for principal windbreak 
(m), Ho is mature height of windbreak (m), and c is a con¬ 
stant, which is determined by the microclimate at the wind¬ 
break location and the type of windbreak. 

However, few studies have quantified the magnitude of 
the constant c and mature height of tree windbreaks (H 0 ) 
because the relationships between these factors and wind 
are very complex (Cao 1983; Jiang 1992). It is usual to 
determine the constant c by considering just wind. Mature 
height is, on the other hand, estimated based on the ex¬ 
periences, or qualitative analysis by the designers. The 
effects of windbreak structure have not been considered in 
equation (1) at all, despite its obvious importance. 

There are several aspects to the problem of determining 
the spacing interval for a system of shelterbelts. The prin¬ 
cipal reason for planting any shelterbelt is to reduce wind 
velocity. The amount and type of wind reduction needed, 
however, depend on the object to be protected and the 
purpose of the protection. Therefore, limiting the investiga¬ 
tion to only one of these purposes does not completely 
solve the problem. For example, Woodruff (1956) found 
that to determine optimum spacing to control soil erosion, 
one must consider the susceptibility to erosion of the soil 
that is to be protected. This varies considerably by soil type 
and makes it difficult to prescribe a general formula for the 
velocity reduction required to control erosion. The object of 
this study was to develop quantitative criteria for the spac¬ 
ing interval of principal tree windbreaks for wind protection 
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based on the mechanical reduction of windspeed, wind¬ 
break height, windbreak structure (optical porosity) and the 
microclimate related to the desired windspeed reduction. 

Materials and methods 

Study area 

The experiments were conducted in the northeast of 
China, including several counties of Changtu, Liaoning 
Province; Nongan, Jiutai, Dehui, Shuangyang and Yushu, 
Jilin Province. The general situation of the shelterbelt site is 
a plain. Soil types are black soil, black meadow soil, black 
calcium soil and black sandy soil (Jiang 1992). The ages of 
shelterbelt range between 10 and 34; the rows of the 
shelterbelts are from 3 to 7. The observations of wind re¬ 
duction for both tree shelterbelts and stubble barriers were 
conducted in Baoli Township and Zhujia Township, Changtu 
County, Liaoning Province. 

Determination of independent variables 

Although windbreak length, shape and width may influ¬ 
ence the shelter effect leeward of windbreaks, tree height 
and windbreak porosity are more important (Cao 1983; 
Heisler and DeWalle 1988; Jiang 1992; Leoffler et al. 1992; 
Zhu et al. 2002). Therefore, mature height and porosity are 
considered as the key windbreak structural indices in de¬ 
termining spacing interval for wind control (wind protection). 
The additional major factor to influence windbreak effec¬ 
tiveness is the problem wind velocity because a reduction 
in momentum implies a reduction in velocity (Plate 1971). 
As described in equation (1), the constant c is mainly de¬ 
termined by microclimate and windbreak type. Thus, we 
defined the spacing interval for principal tree windbreaks 
as: 

SIP = SL rp H o (2) 

Here 8 is the structure index (in this particular study, the 
structure of windbreak is represented by optical porosity), 
Lrp is a parameter of the microclimate, and determined by 
problem winds. The others are the same as mentioned 
above. If we assume 8 Lrp=c, equation (2) becomes equa¬ 
tion (1). The purpose of this particular study is to determine 
quantitatively the three variables in equation (2). 

Windbreaks with various porosities 

The tree windbreaks were composed of poplar (Populus 
spp.), with average windbreak height varied from 8.0 m to 
14.0 m, and the number of rows varied from 3 to 7 rows. 
Optical porosity ( p ) of tree windbreaks with leaves was 
obtained by using the digitizing techniques described by 
Kenney (1987) and Jiang (1989). This method required 
black and white photographic silhouettes of each tree 
windbreak, which was digitized to calculate windbreak op¬ 
tical porosity. The full height porosity was determined by the 
weighted mean of the bottom-half porosity and upper-half 


or crown-half porosity. Windbreak height was defined as 
the distance from the ground surface to the top of the tallest 
tree in the depicted windbreak segment. Nine shelterbelts 
were observed. The full height optical porosities for these 
tree windbreaks were 0.13, 0.15, 0.18, 0,20, 0.23, 0.25, 
0.28, 0.30 and 0.33. 

Reviews of the aerodynamics of barriers by Plate (1971) 
showed that the airflow patterns about wind barriers are 
very complicated. In the simplest case of a barrier normal to 
the airflow, the flow is turbulent (Hagen and Kidmore 1981). 
Barriers are usually within the so-called “constant-stress 
region” near the surface where the open field mean velocity 
profile is described by: 

U mo /U* =(\/K)\ n [{z-D)/Z 0 ]+ f(R,) (3) 

Here Umo is mean velocity at height Z, U- is friction velocity, 
k is von Karman's constant, Z is height of measurement 
from some reference plane, D is zero plane displacement 
height, Z D is a roughness parameter, f(R,) indicates that 
mean velocity also depends on Richardson’s number (R,), 
which is a measure of atmospheric stability. Under neutral 
stability, f(R,) is near 0 (Hagen and Kidmore 1981; Wang 
and Takle 1995, 1996). 

From equation (3), it is obvious that both surface 
roughness and stability of the air influence airflow above 
barriers. In an atmospheric boundary layer with neutral 
stability, the airflow patterns near two geometrically similar 
barriers a and b will be similar in the following way: 

(H/Z 0 ) a =(H/Z 0 ) b (4) 

Here H is windbreak or barrier height (Tanaka et al. 1955; 
Lyles and Allison 1976; Hagen and Kidmore 1981; Zhu 
1981; Murai et al. 1992). This relationship means that 
model barriers that are shorter than the tree windbreaks 
can be used to increase the number of porosities in the 
study. 

Model barriers with known height and optical porosity 
were constructed using one to seven rows of sorghum 
stubbles, which can be re-arranged freely. Porosities of a 
single row were 0.20, 0.40, 0.60 or 0.80. By varying poros¬ 
ity of individual rows and the number of rows, 28 different 
types of model barriers were produced, but because three 
combinations had similar porosity, only 25 kinds of poros¬ 
ities were obtained (Table 1). The maximum height of the 
stubble barriers was 1.85 m tall, and each barrier was 50 m 
long. The model barriers were established in the field, on a 
site that is flat in every direction for at least 2 km. The op¬ 
tical porosity of the stubble barriers was determined by a 
method of dotted grid, which was developed by George et 
al. (1963). 

Windspeed measurements 

Windspeeds and wind directions were measured by us¬ 
ing 3-cup anemometers (528-MODEL and DEM-6 MODEL). 
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The instruments have a starting windspeed of 0.5 m-s' 1 
and 0.2 m-s’ 1 , respectively, and a maximum windspeed of 
30 m-s’ 1 . One-minute average windspeed and direction 
can be read directly from the instruments. The errors of the 
instruments are ±0.2 m-s' 1 and ±0.5 m-s' 1 for DEM-6 
MODEL and 528-MODEL, respectively. Both instruments 
are designed by Meteorological-oceanic Instrument Com¬ 
pany, Tianjin, China. The arrangements for both tree wind¬ 
breaks (DEM-6 MODEL) and the stubble barriers 
(528-MODEL) are as follows: 


Table 1. Average optical porosity with various rows of model 
barriers* 


Porosity of a 




Rows 




single row 

1 

2 

3 

4 

5 

6 

7 

0.20 

0.20 

0.15 

0.10 

0.06 

0.03 

0.01 

0.00 

0.40 

0.40 

0.26 

0.18 

0.13 

0.09 

0.06 

0.04 

0.60 

0.60 

0.47 

0.37 

0.28 

0.21 

0.15 

0.12 

0.80 

0.80 

0.69 

0.61 

0.54 

0.48 

0.43 

0.40 


’Example: In the case of the porosity of a single row is 0.20; the com¬ 
position of model barrier can be formed by 1 row. 2 rows, 3 rows, 4 
rows, 5 rows, 6 rows and 7 rows. The optical porosities of the model 
barriers with 1 row, 2 rows, 3 rows, 4 rows, 5 rows, 6 rows and 7 rows 
are 0.20, 0.15, 0.10, 0.06, 0.03, 0.01 and 0.00, respectively. Three 
combinations produced the same porosity. 

In the case of both tree shelterbelts and stubble barriers, 
for producing the open field values, one instrument was 
placed on the windward side at a distance not less than 20 
times the heights of the windbreaks (one person all through 
read the wind data there during the observed period). The 
distances are expressed in terms of windbreak or model 
barrier height (H). 

In the case of tree shelterbelt, the experiments were con¬ 
ducted during April and June (full leaf period of poplar 
shelterbelt). Windspeed and wind direction were collected 
at -10, -5, -3, -1, 1, 3, 5, 7, 10, 15, 20, 25 and 30 H at a 
height of 2.0 m above the ground; wind data at each point 
were observed simultaneously (Total 13 persons read the 
windspeed and wind direction at an interval of one minute 
at the same time). Data collection (significant data, i.e. 
wind blowing perpendicular to windbreak, and wind velocity 
more than 5 m-s' 1 ) of all anemometers lasted more than 10 
min long at each point, and three repeats were conducted. 
The one minute average windspeed was in range of 
6.0-8.4 m-s’ 1 in the field. 

In the case of stubble barriers, the experiments were 
conducted during October and December, and windspeed 
and wind direction data were collected at -10, -5, -2, -0.25, 
0.25, 5, 10, 15, 20, 25, 30 H from the stubble barriers at 0.5 
m above the ground. Observations were also conducted at 
interval of one min, significant data collection lasted 5 min 
at each point, and three repeats were conducted. The 
one-minute average windspeed was during 5.0-13.0 m-s' 1 
in the field. Wind data at each point were not simultane¬ 
ously observed because of lack of enough instruments. 


Measurements were made only with winds blowing per¬ 
pendicular to the tree windbreaks or the stubble barriers, 
and velocities over 5 m-s' 1 . Average values of 5 min and 10 
min for stubble barrier and tree windbreak respectively at 
per setting were used in the analysis. 

Mature height for tree windbreak 

Measurements on the growth of 118 field tree windbreaks 
were used to determine the mature height for tree wind¬ 
breaks. Age, diameter (DBH) and height of four poplar 
species, Populus xiaozuanrica, P. canadensis, P. simonii 
and P. pseudo-simonii , were measured in the northeast 
China. Stem-analysis techniques were used in an 
age-height growth model for each species. About 80% of 
the tree shelterbelts are Populus xiaozuanrica (69) and P. 
canadensis (24), and 20% are P. simonii (14) and P. 
pseudo-simonii ( 11). 

Both theory, with few exceptions, and empirical meas¬ 
urements have shown that the protected zone associated 
with a windbreak is directly proportional to the height of the 
windbreak (Tanaka et al. 1955; Caborn 1965; Cao 1983; 
Dronen 1988; Finch 1988; Ticknor 1988; Jiang 1994a). The 
height of non-living windbreaks is fixed with their materials 
by designers, but the height of tree windbreak changes with 
the growth of trees. Windbreaks are expected to provide 
permanent and complete wind protection, but it is known 
that the complete protection is not available for several, 
even several tens of years from the planting of windbreaks 
until they reach sufficient height to provide the complete 
protection. 

Some researchers estimate the mature height of tree 
windbreak as the dominant height of mature windbreaks in 
best-adapted species for a given site qualitatively (Cao etal. 
1981, 1983; Wight 1988; Jiang 1992). However, if the es¬ 
timation of mature height is high enough, the spacing in¬ 
terval is too long for some tree species to provide the com¬ 
plete protection (Finch 1988). On the other hand, if the 
mature height is underestimated, the spacing interval for 
windbreaks is short enough, for the complete protection to 
be easily obtained at an early age. Complete protection will 
then quickly extend to the next windbreaks. Therefore, the 
mature height ( H 0 ) is one of the key indices to obtain the 
optimum spacing interval. In this connection, “optimum” 
means smallest proportion of land devoted to tree wind¬ 
breaks and the largest complete protection for the remain¬ 
ing land. The mature heights of tree windbreaks were de¬ 
termined from the investigations and growth development, 
that is, stem analysis techniques. 

There is much information about tree height growth 
models, and the quantitative analysis of forest growth has 
progressed rapidly both in extent and in degree of sophis¬ 
tication (Jerome 1995). Previous studies about the growth 
of trees in shelterbelts (Jiang 1992; Zhu et al. 1992, 1996) 
indicated that the trees in windbreaks have their charac¬ 
teristics in growth because of the special site and edge 
effects for multiple-row windbreaks. Tree growth in most 
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shelterbelts can be well described using the following 
growth model: 

W,(0 = ff M ,{l + P,exp[-y l (f-« b ))r 1 (5) 

where H,{t) is total growth of windbreak height (m), / repre¬ 
sents tree species composed of the windbreaks, t is age of 
windbreak stand (year), f 0 is age at beginning (year), H ias is 
asymptotic growth for a given tree species on a given site 
(m), n, is coefficient related to growth at a particular site and 
y is internal incremental rate for a given tree species. The 
parameters in equation (5) were determined according to 
stem analysis. 

Although relative growth rate and average growth rate 
are usually major indictors for division of tree growth stages, 
the incremental acceleration (A(/)) in growth rate is another 
important index for the growth stage division, particularly for 
the height growth (Husch et al. 1972). The incremental 
acceleration in growth rate in equation (5) is written as: 

A,(t) = H l 2 (t)/to 2 = H m w l \xtf-Y l (t-t 0 )]{ii l 

exrf-y, (t -r 0 )] -1} (1 + Aexrt-K, (/ - fo)] }~ 3 

In equation (6), A(f) has the minimum value in the fol¬ 
lowing format. 

A»j min ~ [lri(Ai) c]/y f (7) 

We found that the tree height of windbreaks would be 
relatively “stable” at age A,™,, that is, the further growth of 
H,{f) after the age of A m „ is relatively slow, thus, the effec¬ 
tive protective distance is “relative stable” because the 
effective protective distance is direct proportion to tf 7 (Ta¬ 
naka et al. 1955). Therefore, we defined the height of tree 
windbreak at age of A m „ as the mature height (H 0 ). 


Results 

Windspeed reduction with various porosities in tree 
windbreaks 

Average windspeeds (three 10-minute averages) at each 
point are listed in Table 2. In this analysis only leeward data 
were used. The reduction of windspeed expressed by 
1 " UieJUopen (Ulee, Uopen are windspeeds in leeward and 
open field respectively) in 0 to 10 H, 0 to 20 H and 0 to 30 H 
of tree windbreaks is shown in Fig.1. From Fig.1, the opti¬ 
mum reduction of windspeed was observed when tree 
windbreak porosity was around 0.25 for all distances lee¬ 
ward of the windbreak. However no significant difference 
between porosity 0.25 and the other various porosities can 
be found at level 0.01 by F-test. 



Fig. 1 Average windspeed reduction of tree windbreaks in 
leeward of 0 to 10 H, 0 to 20 H and 0 to 30 with different po¬ 
rosities, 

(H is height of windbreak or model barrier, and some of the data were 
published in 1992, Jiang 1992) 


Table 2. Average windspeed (Three 10-minute average) at the observed points of tree windbreaks with various porosities (nvs' 1 )* 


Porosity 

Row 

-10H 

-5H 

-3H 

-1H 

1H 

3H 

5H 

7H 

10H 

15H 

20 H 

25H 

30H 

Field 

0.13 

7 

8.14 

7.76 

7.65 

7.35 

7.57 

7.42 

6.49 

5.53 

4.99 

4.81 

6.24 

7.06 

7.61 

7.98 

0.15 

4 

7.69 

7.11 

6.83 

6.52 

7.27 

7.25 

6.03 

5.28 

4.60 

4.40 

5.71 

6.53 

7.21 

7.54 

0.18 

7 

8.62 

7.96 

7.64 

7.29 

8.15 

8.12 

6.73 

5.88 

5.11 

4.89 

6.38 

7.30 

8.07 

8.56 

0.20 

4 

8.53 

7.92 

7.69 

7.02 

8.00 

7.82 

6.37 

6.07 

4.98 

4.50 

5.00 

6.98 

7.88 

8.38 

0.23 

4 

5.95 

5.45 

5.21 

4.95 

5.56 

5.55 

4.50 

4.04 

3.48 

3.29 

4.22 

5.02 

5.55 

6.02 

0.25 

4 

7.07 

6.15 

5.81 

5.32 

6.42 

6.71 

5.06 

4.89 

3.98 

3.87 

4.88 

5.94 

6.66 

7.07 

0.28 

5 

7.38 

6.65 

6.06 

6.13 

6.95 

6.95 

5.42 

4.96 

4.51 

4.17 

5.53 

6.39 

6.89 

7.40 

0.30 

3 

5.75 

5.29 

5.06 

4.82 

5.39 

5.38 

4.41 

3.97 

3.45 

3.27 

4.14 

4.89 

5.38 

5.62 

0.33 

5 


6.75 

6.63 

6.44 

7.05 

6.92 

6.36 

4.60 

4.16 

4.22 

6.35 

5.88 

6.62 

6.99 


*-10H, -5H, -3H, -1H: windward side of the shelterbelt; 1H, 3H, 5H, 7H, 10H, 15H, 20H, 25H, 30H: leeward side of the shelterbelt. 


Windspeed reduction with various porosities in stubble 
barriers 

Horizontal velocity measurements were made in each 
combination of porosity stubble barriers. Because the wind 
data at each point were not measured at the same time, the 
relative windspeed (LWC/ 0 pen) was used in this analysis. 
Mean ratios of leeward to windward windspeeds were 


computed from the cup anemometer at 0 to 5 H, 0 to 10 H, 
0 to 15 H, 0 to 20 H, 0 to 25 H, 0 to 30 H and 0 to 35 H. The 
largest reduction of windspeed was observed in the stubble 
barriers with porosity of 0.13 (Table 3, Fig. 2). The stubble 
barriers with porosities in the range of 0.06 to 0.28 can 
provide better sheltering than the barriers with other poros¬ 
ities. 
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Table 3. Windspeed reduction by stubble barriers with various porosities ( U\JU orxn ) 


No. 

Porosity of No. of 
single line row 

Porosity 

Upwind 

10 H 

Upwind 

5 H 

Leeward 

5 H 

Leeward 
10 H 

Leeward 
15 H 

Leeward 
20 H 

Leeward 
25 H 

Leeward 
30 H 

Leeward 
35 H 

1 

0.20 

1 

0.20 

0.84 

0.74 

0.29 

0.36 

0.46 

0.53 

0.60 

0.64 

0.68 

2 

0.20 

2 

0.15 

0.84 

0,76 

0.18 

0.27 

0.37 

0.45 

0.52 

0.58 

0.63 

3 

0.20 

3 

0.10 

0.84 

0.76 

0.20 

0.31 

0.41 

0.48 

0.55 

0.59 

0.65 

4 

0.20 

4 

0.06 

0.80 

0.71 

0.19 

0.31 

0.39 

0.47 

0.53 

0.58 

0.63 

5 

0.20 

5 

0.03 

0.81 

0.72 

0.21 

0.32 

0.41 

0.48 

0.54 

0.60 

0.65 

6 

0.20 

6 

0.01 

0.84 

0.75 

0.26 

0.36 

0.45 

0.52 

0.58 

0.65 

0.69 

7 

0.20 

7 

0.00 

0.86 

0.77 

0.28 

0.38 

0.46 

0.54 

0.61 

0.66 

0.72 

8 

0.40 

1 

0.40 

0.92 

0.88 

0.51 

0.54 

0.60 

0.65 

0.70 

0.73 

0.77 

9 

0.40 

2 

0.26 

0.85 

0.79 

0.25 

0.31 

0.40 

0.48 

0.54 

0.59 

0.63 

10 

0.40 

3 

0.18 

0.86 

0.80 

0.16 

0.27 

0.37 

0.45 

0.52 

0.58 

0.63 

11 

0.40 

4 

0.13 

0.82 

0.76 

0.13 

0.24 

0.35 

0.43 

0.50 

0.55 

0.60 

12 

0.40 

5 

0.09 

0.83 

0.76 

0.17 

0.29 

0.38 

0.46 

0.53 

0.58 

0.63 

13 

0.40 

6 

0.06 

0.83 

0.76 

0.17 

0.28 

0.37 

0.45 

0.52 

0.57 

0.62 

14 

0.40 

7 

0.04 

0.86 

0.79 

0.22 

0.34 

0.42 

0.49 

0.54 

0.59 

0.64 

15 

0.60 

1 

0.60 

0.92 

0.89 

0.58 

0.65 

0.71 

0.76 

0.80 

0.83 

0.85 

16 

0.60 

2 

0.47 

0.90 

0.87 

0.41 

0.47 

0.54 

0.60 

0.65 

0.68 

0.72 

17 

0.60 

3 

0.37 

0.87 

0.83 

0.32 

0.39 

0.47 

0.54 

0.59 

0.63 

0.67 

18 

0.60 

4 

0.28 

0.88 

0.84 

0.25 

0.34 

0.43 

0.50 

0.57 

0.62 

0.66 

19 

0.60 

5 

0.21 

0.84 

0.78 

0.23 

0.33 

0.42 

0.49 

0.57 

0.60 

0.64 

20 

0.60 

6 

0.15 

0.88 

0.81 

0.22 

0.31 

0.41 

0.49 

0.54 

0.52 

0.62 

21 

0.60 

7 

0.12 

0.87 

0.82 

0.23 

0.35 

0.44 

0.51 

0.57 

0.62 

0.67 

22 

0.80 

1 

0.80 

0.89 

0.88 

0.82 

0.83 

0.85 

0.87 

0.89 

0.90 

0.91 

23 

0.80 

2 

0.69 

0.87 

0.86 

0.71 

0.73 

0.76 

0.79 

0.81 

0.82 

0.84 

24 

0.80 

3 

0.61 

0.89 

0.86 

0.71 

0.71 

0.74 

0.77 

0.79 

0.82 

0.84 

25 

0.80 

4 

0.54 

0.88 

0.86 

0.66 

0.68 
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26 

0.80 

5 

0.48 
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0.59 
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0.78 

27 

0.80 

6 
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0.84 
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28 

0.80 

7 

0.40 

0.90 

0.87 

0.53 

0.57 
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0.67 

0.71 

0.74 

0.77 



Windbreak porosity 


Fig. 2 Mean ratios of leeward to windward windspeeds of 
windbreak barriers with various porosities in different range 
of locations. 


Mature height 

The number of trees used for stem-analysis: Populus 
xiaozuanrica, 15; P. canadensis, 6; P. simonii, 5 and P. 
pseudo-simonii , 8, which were used to develop the 
age-height growth model in equation (5). The data from the 
model were compared with the field measurements, and no 
significant differences were found. The results of this 
analysis have been used in Changchun District, Jilin Prov¬ 
ince, China, and their reliability has been proven in a pre¬ 
vious study (Jiang, 1996). The mature heights of four poplar 
species were determined according to equations (5), (6), 


and (7) and listed in Table 4. 

Table 4. Determination of mature height (H 0 ) of poplar wind¬ 
breaks 


Growth parameters A vmn Mature height 


Tree species 

Hias A* i 

Y< 

fl 2 

/a 

(Ho) /m 

Populus simonii 

17.9 

8.16 

0.159 

0.828 

23.5 

13.4 

P. pseudo-simonii 

19.9 

9.40 

0.154 

0.901 

24.0 

15.6 

P. canadensis 

25.7 

5.80 

0.206 

0.867 

15.0 

20.3 

P. xiaozuanrica 

21.0 

13.77 

0.239 

0.839 

15.5 

17.5 


Generalized growth model is equation (4), A vmn as de¬ 
fined in equation (6), F? is coefficient of determination, and 
each model was significant at p <0.01 

Discussion 

Structure index S 

Wind reduction is closely related to windbreak porosity, 
and classifying tree windbreaks on the basis of porosity 
would facilitate comparisons of effectiveness among tree 
windbreaks (Hagen and Kidmore 1971; Jiang etal. 1994b). 
The optimum aerodynamic porosity for field shelterbelts is 
usually considered to be 0.35 to 0.45 (Hagen and Kidmore 
1971; Loeffler etal. 1992), but unfortunately, it is impossible 
to physically measure the aerodynamic porosity of field tree 
windbreaks. Much effort has therefore been directed to¬ 
wards finding alternative measurements. Optical porosity 
( 0), a two-dimension measure of porosity determined from 
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the windbreak silhouette, may be a promising alternative to 
aerodynamic porosity. Heisier and DeWalle (1988) suggest 
that the optical porosity is a useful guide to windbreak 
functions, and therefore an appropriate guide for evaluating 
narrow windbreaks. 

The optimum optical porosity has been studied by many 
researchers (Tanaka et al. 1955; George et al. 1963; Moy¬ 
sey and Mcpherson 1966; Hagen and Kidmore 1971,1981; 
Zhou et al. 1991; Jiang 1992), but there is considerable 
variation in results. The source of this variation includes 
differences in shelterbelts structure, the effects of thermal 
instability in the field, the type of instruments used, the tree 
species, windbreak width, and the method used to deter¬ 
mine the optical porosity. Despite this variation, most stud¬ 
ies have suggested that the best sheltering is obtained 
using windbreaks with porosities of 0.20 to 0.50 (Moysey 
and Mcpherson 1966; Dronen 1988; Jiang 1992; Jiang etal. 
1999). 

The results in this experiment indicated that the optimum 
porosities measured in the tree windbreaks and in the 
stubble windbreak barriers were not consistent. Although 
similarity was retained between the two experiments using 
tree windbreaks and stubble barriers, we recognize that the 
stubble barriers do not react exactly like the tree wind¬ 
breaks. Generally, under the tree windbreak conditions, the 
branches, twigs and leaves of trees move with the wind, i.e., 
the porosities of tree windbreaks are changing continuously 
with the wind velocity. However, the porosities of the stub¬ 
ble barriers seldom change, except with extremely strong 
winds. The tendencies of windspeed changes with different 
porosities around the optimum values (around 0.13 for the 
barriers, and 0.25 for the tree windbreaks) are very similar 
between the tree windbreaks and the stubble barriers (Fig.1 
and Fig.2). These results, both from the tree windbreaks 
and the stubble barriers, are very similar to Moysey and 
Mcpherson’s (1966). Those authors suggested that wind¬ 
breaks with porosities in the range of 0.15 to 0.30 have 
better sheltering, and approximately 0.25 porosity is rec¬ 
ommended for the best sheltering. Therefore, The regres¬ 
sion relationship between windbreak porosity and wind- 
speed reduction, namely, efficiency of windbreak in wind 
reduction changing with porosity, which was defined as 
structure index 8 , can be obtained based on the observa¬ 
tions of wind reduction from tree shelterbelts with combin¬ 
ing the results from stubble barriers (Fig. 3). The multiple 
regression equation was: 

8 =5.04249 p 3 -8.7712 P 2 +3.4239 /3 +0.6139 

(0< P <0.80) (8) 

Ff =0.9339, significant at p<0.01 

Here P is windbreak porosity, it varies from 0.00 to 0.80, 
the range of value S is from 0.31 to 1.00. 

When p =0.00, S =0.61, when p =0.80, <5 =0.31 and 
when p =0.25, 8 =1.00. 



Porosity 

Fig. 3 Efficiency of windspeed reduction for field windbreaks 
with various porosities 

(The best wind reduction is observed from a windbreak with porosity of 0 25) 

Parameter of problem wind L rp 

The most frequent design requirement is that the wind be 
reduced below the dangerous level over a maximum dis¬ 
tance. The reduction in wind velocity needed to protect 
objects from damage depends on the climatology of the 
approach wind conditions such as velocity of open wind, 
direction of problem wind and temperature, soil type and 
water content, the methods of tillage used and the situation 
of protected objects. The literature indicates that threshold 
velocity (threshold velocity defined as: velocity necessary to 
start movement of a soil particle or initially to damage the 
protected objects) is a key index for determining the effec¬ 
tiveness of windbreaks (George etal. 1963; Cao etal. 1981; 
Dronen 1988; Finch 1988; Wight 1988). 

If problem wind velocity is V pv/ , and the threshold velocity 
is Vm (V pw >Vih), the reduction of windspeed U m can be 
written as: 

U m = (V -V,*)/V100% = 100 - V th IV pw \ 00% (9) 

Based on the results of wind reduction under various 
porosities, the horizontal wind-profile with optimum porosity 
can be obtained as shown in Fig.4. The wind velocity re¬ 
duction with the windbreak under the optimum porosity 
(0.25) from 5 to 35 H expressed as percentage of the open 
wind can be well described by a logarithm form as: 

V lee /U open =V' n ( L r P ) + <l>Q ( 10 ) 

where V, o 0 are empirical constants, the others are the 
same as mentioned above. 

In equation (10), Lrp is the only independent variable 
that significantly contributed to predicting the relative 
windspeed ( Uie</U op en ), namely, the reduction of windspeed 
of problem wind ( U m ) is determined by L v . Therefore, inte¬ 
grating or combining equations (9) and (10), i.e., replacing 
Vth/V P w with Uiee/Uo P en, the relationship regression is shown 
in equation (11). 

U m =m-yln{L rp ) + (i> 0 (11) 

Solving equation (11) with the data from Fig.4, Lrp can be 
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written as: 

L rp = exp[(130.3338-(/ w )/34.2176J (12) 

fl 2 =0.9949, significant at pcO.OI 



Relative protected distance (H) 

Fig. 4 Standard horizontal windspeed-profile of windbreaks 
with optical porosity of 0.25 

It should be noted that the determination of the threshold 
velocities ( Vm) for various protected objects is very complex, 
and equations (10) and (11) were obtained from the range 
of 5H to 35H of windbreaks or stubble barriers, therefore, 
the value of U m is limited between 9% and 75%. The range 
of Um is determined based on the calculation by equation 
(12). The result is only limited to wind protection, not in¬ 
cluding snow accumulation and other purposes. 

Summary and conclusions 

Field investigations were conducted to obtain information 
on the most important factors for determining the spacing 
interval. The results suggest that the average minimum 
wind reduction occurred in stubble barriers with 0.13 po¬ 
rosity (0 to 35 H), and tree windbreaks with 0.25 porosity (0 


to 30 H). The stubble barriers with porosities in the range of 
0.06 to 0.28 provide better sheltering than those with other 
porosities. The degree of efficiency of windbreaks depends 
not only on windbreak structure (optical porosity), but also 
on the type of protection desired or the problem winds. The 
influence of shelterbelt structure (optical porosity) upon 
leeward windspeed reduction is given by a multiple regres¬ 
sion equation based on the observed results. Leeward 
relative windspeed reduction as an important sheltering 
characteristic was used to determine the parameters of 
structure index (<5 ) and problem wind (L v ) for the spacing 
interval. The mature height of tree windbreak was also 
determined by developing a growth model of trees in wind¬ 
breaks. 

As a prerequisite to evaluating and comparing the spac¬ 
ing interval of the different porosity combinations of 
shelterbelts, it is necessary to establish desired levels of 
velocity reduction. As an example, levels of reduction of 
relative windspeed (Urw) varying from 10% to 70% were 
chosen to calculate parameter of problem wind (Lrp). While 
these are only arbitrary levels based on the result of this 
particular study, one or the other should be appropriate for 
most wind velocity reduction problems. Porosities of wind¬ 
breaks varying in some ranges of 0.1 to 0.45, 0.12 to 0.40, 
0.15 to 0.36 and the optimum of 0.25 were chosen to cal¬ 
culate the structure index (8). Mature heights ( H 0 ) of four 
poplar species were given in Table 4 (Poplar shelterbelts). 
The results showing spacing intervals for the principal 
poplar windbreak are given in Table 5. All the results dis¬ 
cussed were based on the measurements of windspeeds in 
plain farmland, the relationship examined are therefore not 
applicable to other disaster protection such as snow man¬ 
agement. 


Table 5. Spacing interval (SIP) for principal poplar tree windbreaks (m)* 


Tree species 

Range of porosity 

(P) 

Average of 
Index (5) 


Reduction of relative windspeeds required 


UrwIO 

Urw20 

Urw30 

Urw40 

Upw50 

U rw60 

Urw70 

Populus simonu 

0.10 to 0.45 

0.948 

426.2 

318.2 

237.5 

177.3 

132.4 

98.8 

73.8 


0.12 to 0.40 

0.966 

434.3 

324.2 

242.1 

180.7 

134.9 


75.2 


0.15 to 0.36 

0.981 

441.0 

329.2 

245.8 

183.5 

137.0 

102.3 

76.4 


Optimum 0.25 

1.000 

449.5 

335.6 

250.6 

187.1 

139.7 

104.3 

77.8 

P.pseudosimonii 

0.10 to 0.45 

0.948 

497.4 

371.3 

277.2 

207.0 

154.5 

115.4 

86.1 


0.12 to 0.40 

0.966 


378.4 

282.5 

210.9 

157.5 

117.6 

87.8 


0.15 to 0.36 

0.981 

514.7 

384.2 

286.9 

214.2 

159.9 

119.4 

89.1 


Optimum 0.25 

1.000 

524.6 

391.7 

292.4 

218.3 

163.0 

121.7 

90.9 

P. canadensis 

0.10 to 0.45 

0.948 

646.8 

482.9 

360.5 

269.1 

200.9 

150.0 

112.0 


0.12 to 0.40 

0.966 


492.0 

367.3 

274.2 

204.7 

152.9 

114.1 


0.15 to 0.36 

0.981 

669.3 

499.7 

373.0 

278.5 

207.9 

155.2 

115.9 


Optimum 0.25 

1.000 

682.2 

509.3 

380.3 

283.9 

212.0 

158.2 

118.1 

P. xiaozuanrica 

0.10 to 0.45 

0.948 

558.6 

417.1 

311.4 

323.5 

173.6 

129.6 

96.7 


0.12 to 0.40 

0.966 

569.3 

425.0 

317.3 

236.9 

176.9 

132.0 

98.6 


0.15 to 0.36 

0.981 

578.1 

431.6 

322.2 

240.6 

179.6 

134.1 

100.1 


Optimum 0.25 

1.000 

589.3 

440.0 

328.5 

245.2 

183.1 

136.7 

102.0 


’Spacing interval (m) is from equation (2), 8 defined as equation (8), LUto to U ra7 o are the values of reduction of relative windspeeds required, U m , 0 


to Uw7o represent 10% to 70% reduction of relative windspeeds, L v is obtained from equation (12) by using U, w )0 to U„?o', mature height (H 0 ) is from 
Table 2. 


Considerable further study is needed to determine the reported here is indicative that spacing interval ( SIP) is 

threshold velocities for various protected objects. The study directly related to relative windspeed reduction (U m ), but 
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not the threshold velocity. 

Acknowledgement 

We'd like to thank the people who attended the research 
programs presided over by Prof. Feng-qi Jiang during 1986 
and 1996 for their great help in material preparations. We 
gratefully acknowledge Dr. Roest Naomi M. for his valuable 
comments and suggestions on the manuscript. Thanks are 
also due to the editors and referees for their helpful com¬ 
ments and careful revision on this manuscript. 

References 

Caborn, J.M. 1965. Shelterbelts and microclimate [Mj. Faber and 
Faber ltd., London, p48-169. 

Cao Xinsun, Lei Qidi, Jiang Fengqi. 1981. Optimum porosity and 
transversal section for shelterbelts [J], Bull. Inst. For. Pedo., 
Academia Sinica, 5: 9-19. (In Chinese) 

Cao Xinsun. 1983. Shelterbelts for farmland [M]. Beijing: Chinese 
Forestry Press, p98-189. (In Chinese) 

Dronen, S.l. 1988. Layout and design criteria for livestock wind¬ 
breaks [J]. Agri., Ecosys. Envir., 22/23: 231-240. 

Finch, S.J. 1988. Field windbreaks: design criteria [J]. Agri., Eco¬ 
sys. Envir., 22/23: 215-228. 

George, E.J., Broberg, D. and Worthington, E.L. 1963. Influence of 
various types of field windbreaks on reducing wind velocities and 
depositing snow [J], J. For., 61: 345-349. 

Hagen, L.J. and Kidmore, E.L. 1971. Windbreak drag as influenced 
by porosity [J], Trans. ASAE, 14: 464-465. 

Hagen, L.J. and Kidmore, E.L. 1981. Simulation of effect of wind 
barriers on airflow [J]. Trans. ASAE, 24: 1002-1008. 

Heisler, G.M. and DeWalle, D.R. 1988. Effects of windbreak struc¬ 
ture on wind flow [J]. Agri., Ecosys. Envir., 22/23: 41-69. 

Husch, 8., Miller, C.l. and Beers, T.W. 1972. Forest measuration. 

The Ronald Press Co.. New York, p281 -319. 

Jerome, K.V. 1995. Synthesis- growth models for tropical forests, a 
synthesis of models and method [J], For. Sci., 41: 417-42. 

Jiang Fengqi, Fu Menghua, Xu Jiyan. 1989. Estimating windbreak 
porosity using digitized photographic silhouettes [C]. In: Xiang 
Kaifu el al. (eds) Protective plantation in Northeast of China. 
Harbin: Northeast Forestry University Press, p399-401. (In Chi¬ 
nese) 

Jiang Fengqi. 1992. Techniques and theoretical foundation of 
management for shelterbelt. Beijing: Chinese Forestry Press, 
pi-29. (In Chinese) 

Jiang Fengqi, Zhou Xinhua, Fu Menghua, Zhu Jiaojun. 1994a. 
Shelterbelt porosity model and its application [J], Chin. J. Appli. 
EcoL, 5: 251-255. (In Chinese) 

Jiang Fengqi, Zhu Jiaojun, Zhou Xinhua, Lin Heming. 1994b. Pro¬ 
tective maturity and regeneration of shelterbelt [J]. Chin. J. Appli. 
Ecol., 5: 337-341. (In Chinese) 

Jiang Fengqi. 1996. Rational management and improving tech¬ 
niques for protective forests. Beijing: Chinese Forestry Press, 
pi-6. (In Chinese) 

Jiang Fengqi, Zhu Jiaojun, Zhou Xinhua. 1999 Model of continu¬ 
ous economic effects of shelterbelts or windbreaks and its 
applications [J], Scientia Silvae Sinicae, 35: 16-21. 

Kenney, W.A. 1987. A method for estimating windbreak porosity 
using digitized photographic silhouettes [J]. Agri. For. Meteo., 39: 
91-94. 


Loeffler A.E., Gordon, A.M, and Gillespie, T.J. 1992. Optical po¬ 
rosity and windspeed reduction by coniferous windbreaks in 
Southern Ontario [J]. Agrofor. Sys., 17: 119-133. 

Lyles, L and Allison, E. 1976. Wind erosion: The protection role of 
simulated standing stubble [J], Trans, the ASAE, 19: 61-64. 

Maki, T. 1985. Studies on the windbreak nets (9), Variations of 
turbulent characteristics by two and three successive windbreak 
nets [J], J. Agri. Meteo., Jpn, 41: 17-24. 

Moysey, E.B. and Mcpherson, F.B. 1966. Effect of porosity on 
performance of windbreaks [J]. Trans. ASAE, 9: 74-76. 

Murai, H., Ishikawa, M., Endo, J. and Tadaki, R. 1992. The coastal 
forest in Japan. Soft Science, INC., Tokyo, p315-408. (In Japa¬ 
nese) 

Plate, E.J. 1971. The aerodymics of shelter belts [J]. Agri. Meteo., 
8: 203-222. 

Skidmore, E.L. and Hagen, L.J. 1977. Reducing wind erosion with 
barriers [J], Trans. ASAE, 20: 911-915. 

Tanaka, S„ Tanizawa. T. and Sano, H. 1955. Studies on the wind in 
front and back of the shelter-hedges (5) [J], J. Agri. Meteor., Jpn., 
11: 91-94. 

Tibke, G. 1988. Basic principles of wind erosion control [J]. Agri., 
Ecosys. Envir., 22/23: 103-122. 

Ticknor, K.A. 1988. Design and use of field windbreaks in wind 
erosion control system [J], Agri., Ecosys. Envir., 22/23:123-132. 

Wang, H. and Takle, E.S. 1995. A numerical simulation of bound¬ 
ary-layer flows near shelterbelts [J], Boundary-Layer Meteor., 75: 
141-173. 

Wang, H. and Takle, E.S. 1996. On three-dimensionality of 
shelterbelt structure and its influence shelter effects [J], Bound¬ 
ary-Layer Meteor., 79: 83-105. 

Wight, B. 1988. Farmstead windbreaks [J]. Agri., Ecosys. Envir., 
22/23: 261-280. 

Woodruff, N.P. 1956. The spacing interval for supplemental 
shelterbelt [J], J. For., 54: 115-122. 

Zhou Xinhua, Jiang Fengqi, Zhu Jiaojun. 1991. Study on random 
error of shelterbelt porosity estimate by measuring photo with the 
help of digitized photographic silhouettes [J], Chin. J. Appli. Ecol., 
2:193-202. (In Chinese) 

Zhu Jiaojun, Jiang Fengqi. 1992. Management basics and sus¬ 
tainable regeneration pattern for shelterbelts or windbreaks [C], 
In: Proceeding of China association for science and technology, 
First annual meeting of youth, volume agriculture. Beijing: Sci¬ 
ence and technology press of China, 76-83. (In Chinese) 

Zhu Jiaojun, Jiang Fengqi, Zhou Xinhua, Lin Heming. 1996. 
Growth stage and classification of poplar windbreaks [J]. Chin. J. 
Appli. Ecol., 7: 11-14. (In Chinese) 

Zhu, J.J., Matsuzaki, T., Sakioka, K. and Yamamoto, M. 1998. 
Windspeed in a coastal forest belt of Japanese black pine 
—Vertical wind profile [J]. Trans. Jpn. For. Soci., 109: 421-424. 

Zhu, J.J., Matsuzaki, T. and Sakioka, K. 2000. Wind speeds within 
a single crown of Japanese black pine (Pinus thunbergii Pari.) [J], 
For. Ecol. Manag., 135: 19-31. 

Zhu, J.J., Matsuzaki, T. and Gonda, Y. 2001. Wind profiles in a 
coastal forest of Japanese black pine (Pinus thunbergii Pari.) 
with different thinning intensities [J], J. For. Res., the Jpn. For. 
Soc., 6: 287-296. 

Zhu, J.J., Gonda, Y., Matsuzaki, T. and Yamamoto, M. 2002. Salt 
distribution in response to optical stratification porosity and rela¬ 
tive windspeed in a coastal forest. Agrofor. Sys., in press. 

Zhu Tingyao. 1981. Research on wind protection of windbreak in a 
wind tunnel [J]. Bulle. Insti. For. Pedo!., Academia Sinica, 5: 
29-45. 



